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Рассматриваются возможности применения современных пакетов блочно-модульного моделирования в зада- 
чах имитации переходных процессов в газовом потоке на участке транспортировки природного газа. Предла- 
гается подход к моделированию в векторном формате пространственно распределённого объекта ячеечной 
моделью, построенной в специальном нелинейном базисе. Это позволяет решать задачи диагностики и 
управления потоком в переходных процессах за счёт макроячеечного моделирования баланса его массы и 
его кинетической энергии. Адекватность начального и конечного состояний ячеечной модели обеспечивает- 
ся расчётом вектора опорных параметров, характеризующих стационарное течение потока, и сравнением его 
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Введение. Для эффективного управления магистральными потоками природного газа, а также 
диагностики и прогнозирования состояния транспортной системы особое значение имеет возмож- 
ность математического или имитационного моделирования объекта. Это, прежде всего, касается 
задачи моделирования динамических процессов транспортировки газа в реальном времени [1]. 
Процесс транспортировки на участке магистрального газопровода описывается характеристиками 
потока газа, который в связи с пространственной протяжённостью по свойствам можно отнести к 
сплошной среде. Подобные процессы описываются в теории гидромеханики, гидро- и газовой ди- 
намики, а уравнения, описывающие эти процессы, являются уравнениями в частных производных, 
определённых для двумерных функций пространства и времени. 

Анализ публикаций показал, что рассмотренные в монографиях математические модели 
газовой среды либо являются слишком общими, либо приведены к какому-то конкретному виду 
воздействия, например перемещение поршня в сосуде с газом. Поэтому актуальной задачей явля- 
ется получение и исследование применимости математической модели процесса транспортировки 
технологического агрегата, при минимальном количестве допущений; построение алгоритма, до- 
пускающего моделирование реакции на произвольные граничные и параметрические воздействия. 

Серьёзной проблемой компьютерного моделирования являются математические особенно- 
сти модели, такие как нелинейность, пространственно-временная многомерность, большое коли- 
чество дифференциально-интегральных связей, определяющих необходимость применения чис- 
ленных методов приближённого интегрирования. 

Численное моделирование различных технологических процессов всё больше и больше 
применяется для выявления нештатных ситуаций в моделируемом объекте, уменьшения количе- 
ства брака при протекании какого-нибудь технологического процесса. Наибольшую техническую 
завершённость получили те технологии моделирования, которые основаны на моделях математи- 
ческой физики и газогидродинамики, т. е. на тех моделях, которые в докомпьютерный период 
были реализованы как лабораторные установки. Большое количество теоретических исследова- 
ний, в частности изучение свойств моделей при предельных переходах, были проведены и опуб- 
ликованы для моделей сплошной среды — газа, жидкости, кристаллов, электромагнитного поля. 
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Ещё одной развитой областью приложения численного моделирования являются важные для на- 
родного хозяйства технологические процессы. К наиболее законченным решениям в этой области 
относятся компьютерные модели струн, мембран, балок, тепловых полей в твёрдых веществах, 
линий электрической и оптической связи, потоков, обтекающих летательные аппараты и плав- 
средства. 

Моделирование накопительно-транспортировочных переходных процессов, протекающих 
в газотранспортной сети, недостаточно хорошо изучено. Большая часть публикаций посвящена 
моделированию статических режимов, оптимизации режимов компрессорных агрегатов и методам 
диагностики самого газопровода как аппарата транспортировки. В данной работе моделируются 
динамические процессы перераспределения плотности, скорости и давления газа в результате 
переноса и распространения колебаний в движущемся по газопроводу газе. 
Теоретический анализ задачи. Исходная математическая модель, представляющая изотерми- 
ческий процесс (режим, характерный для трубопроводов дальней транспортировки газа при его 
грунтовой прокладке) передачи потоком газа массы вещества и преобразования его механиче- 
ской энергии, описывается математической системой следующего вида: 


д 0 
Е 1 
ОЕ ОХ — 
2 
Е ее Ва ПА Ва (2) 
ОЕ р дх р Н 0х ПО р 








Здесь р = р (х, В, и(х Ви =& (х В =р(х, В и(х 8 — плотность и скорость транс- 
портируемого газа, а также переменная состояния (ПС) потока газа, имеющая физический смысл 
массового расхода, соответственно. Обе переменные изменяются в общем случае в пространстве 
и во времени. Константы модели и начальные условия для переменных состояния заданы сле- 
дующими значениями: А = 8,3144 Дж/моль К; То = 303 К; [= 0,01604 моль; ро =р (0, 0) = 
= 47,7521 кг/м3; = = 506,268 м"; р = 1,22 м; Л = 0,000661 м“/с? [11. 

Рассмотрим подход к дискретизации параметров р = р (х, Ви= = &(х, В модели (1) — (2) 
на примере более простой модели переноса концентрации вещества [2]. Разобьём процесс проте- 
кания на элементарные объёмы — ячейки (см. рис. 1). Допустим, что в пределах ячейки концен- 
трация постоянна, тогда концентрацию на каждом участке можно характеризовать одной пере- 
менной — Си сопоставить её с одной точкой внутреннего пространства участка. Множество таких 
точек называется пространственной сеткой. 

Пространственно-временная аппроксимация градиентов в узлах этой сетки (а иногда и 
между ними), позволяет получить аппроксимированную математическую модель, описываемую 
системой дифференциальных уравнений, которые можно решить с помощью различных числен- 
ных методов. Такой подход, в 
котором сетка строится по 
пространству, а во времени 
модель сохраняет свои конти- 
нуальные свойства, называет- 
ся полудискретизацией. 

Возможности — совре- 
менных пакетов программ в 
плане программной реализа- 
ции алгоритмов, построенных 
на описанных методах, позволяют не только использовать стандартные методы из библиотек 
подпрограмм инженерных расчётов, но и осуществить блочно-модульное построение модели по- 
добного распределённого объекта. Наиболее эффективной и компактной может быть ячеечная 





Рис. 1. Пространственная сетка и расчётная схема моделирования процесса пе- 
редачи концентрации вещества 
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модель любого распределённого процесса в среде пакета УитийпК за счёт использования возмож- 
ностей межблочной передачи векторных сигналов и применения синтаксиса языка системы 
МАТТАВ для дополнительной обработки подобных векторных сигналов. 


го_1 

То \оферасе 
_1 

То \Моферасе1 
ер=_Т 

То \оферасет 

То 'Лоферасе:3 






ерз 


ый 





с сч 

ро р 

с с 

= = ма 
ше ЕТ ше РТ ше а 
|= а СОЯ 2 ма 
|= |= |=: 
5 | 5 Е — 5 |= 
ЕБь от Ешь |от Е от 

-- = 

= = "ша 

и ии 1 

с с 


Чеер=у9Т 


(тои 9+ 


(2 азРуре|пе 
бпапис Мо 9 е | у 


пан ап 9 

Рип ой 

О А 
сепга| 9№егепсе 


арргохитаНоп 
РипеНоп 
СА мВ 


пайА ап9 
арргожитаНоп 


— 
ри 
ре 
= 
р 
— 
ри 

т 

2 

Г. 
— 

= 
= 
ри 
= 





ты с. 
ИЕ — нЕ 


Рис. 2. Модель распределённого объекта в среде пакета ЗититпкК 


СИоск 


Например, при реализации ячеечной модели участка магистрального газопровода (1) — 
(2) элементы этого языка позволяют непосредственно в схеме вычислять конечные разности для 
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аппроксимации градиентов (рис. 2) в (1) — (2). Эти же возможности позволяют обеспечить ини- 
циализацию векторного интегратора одномерным вектором начального состояния. 

В состав схемы моделирования входят: 1 — программный формирователь ступенчатого 
параметрического перехода; 2 — мультиплексор, используемый для сведения значений гранич- 
ных параметров с сигналами о состоянии для вычисления градиентов; 3 — дифференциатор пер- 


вого порядка точности, основанный на следующей формуле вычисления приращения между дис- 
кретными аргументами: 
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Рис. 3. Алгебраические законы формирования изменений на интеграторах, моделирующих состояние 
динамической системы 
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Рис. 4. Переходные процессы по расходу в сечениях участка магистрального газопровода 


Блок 4 реализует моделирование конкретных законов накопления массы, импульса, энер- 
гии (рис. 3) при движении, соотнесённых с узлами сетки, обеспечивая сохранность этих парамет- 
ров в движущемся потоке газа. Мультиплексор 5 из двух векторных сигналов об изменении со- 
стояния формирует один общий, который затем интегрируется блоком 6. Делитель 7 вычисляет 
сигнал, характеризующий скорость потока газа. Группа блоков 8 обеспечивает выборку каждого 
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десятого сигнала для наблюдения на осциллографах 9 и создания истории изменения состояния в 
рабочей области системы МАТИАВ 10. 

На рис. 3 показан модуль схемы, осуществляющий аппроксимацию градиентов, входящих 
в (1), (2), перерасчёт их по этой модели в величины, характеризующие интенсивности накопле- 


ния, и, в соответствие с моделью (1) — (2), равных производным параметров р и 5. 
48 
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Рис. 5. Переходные процессы по плотности в сечениях участка магистрального газопровода 
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Рис. 6. Переходные процессы по скорости потока в сечениях участка магистрального газопровода 


Наиболее исследованными параметрами распределённых объектов являются их статиче- 
ские профили, которые измеряются датчиками и моделируются в контурах АСУ ТП устройствами 
промышленной автоматики. Поскольку моделируемую систему можно рассматривать как сложный 
нелинейный объект, то его статический профиль отождествляется с его аттрактором (притяги- 
вающим многообразием). В силу процессов самовыравнивания плотности и диссипационных про- 
цессов, динамические процессы по переменным, характеризующим состояние объекта, должны 
затухать. Именно такое поведение компьютерной модели свидетельствует о том, что погрешности 
аппроксимации и вычислений не накапливаются, а нивелируются модельной вязкостью и дисси- 


64 


Вестник ДГТУ. 2012. №1 (62), вып. 2 


пацией. Из сказанного можно сделать вывод, что установившееся значение параметров модели 
динамики должно соответствовать конкретизированным уровням статических профилей. 
Результаты исследований. Эксперименты с однонаправленными и разнонаправленными раз- 
ностями, аппроксимирующими градиенты расхода и плотности, показали, что можно добиться 
устойчивого математического моделирования ударных волн и явления массопереноса, если в 
уравнении баланса массы (1) использовать противопотоковые (левосторонние) пространственные 
аппроксимации, а для уравнения импульса разности по потоку — правосторонние. 

Имитационные исследования показали устойчивость переходных процессов, которые схо- 
дятся к уровням, соответствующим профилю, просчитанному после изменения расхода на вход- 
ном участке и плотности на выходном. Графики (рис. 4—6) построены для случая 75 % уменьше- 
ния плотности скачком в момент времени Ёоздействия = 1 С. 

В результате получена компьютерная модель переходных и установившихся режимов 
транспортировки природного газа. Эта модель позволяет провести исследование применимости 
многосеточных методов моделирования, основанных на том свойстве модели, что переходные 
процессы разделяются на две группы — быстрые и медленные (см. рис. 4—6). Данная модель по- 
зволяет создать и протестировать алгоритмы моделирования с повышенным по пространству по- 
рядком точности, за счёт применения методов сплайн-дифференцирования — это позволит ис- 
пользовать более грубую пространственную сетку. 
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